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We present a lattice calculation of the K ^ n and K matrix elements of the AS' — 1 effective weak 
Hamiltonian which can be used to determine e' /e and the AI = 1/2 rule for K decays in the Standard Model. 
The matrix elements for K —> tttv decays are related to A" tt and K using lowest order chiral perturbation 
theory. We also present results for the kaon B parameter, Bk- Our quenched domain wall fermion simulation 
was done at /3 = 6.0 (a~^ ^ 2 GeV), lattice size 16^ x 32 x 16, and domain wall height Ms = 1.8. 



1. Introduction 

Recent measurements of direct CP violation 
(e'/e 7^ 0) in if ^ TTTT decays at FNAL and CERN 
allow an important test of the Standard Model (in 
particular, the CKM mixing paradigm). The ef- 
fective weak Hamiltonian governing strangeness 
changing K decays has been computed to next- 
to-leading order in QCD and QED by the Mu- 
nich and Rome groups; the remaining piece of 
the puzzle is the hadronic matrix elements of the 
operators of this effective weak Hamiltonian. 

The recent advance of domain wall (and over- 
lap) fermions which maintain chiral symmetry 
to a high degree of accuracy allows for a 
new attempt at this old problem. Chiral sym- 
metry of domain wall fermions provides a sig- 
nificant advantage when computing light quark 
QCD observables since the lattice artifacts that 
arise when this symmetry is explicitly broken 
are greatly reduced. Mixing and renormaliza- 
tion of operators, which is already complicated 
in the continuum, is readily handled with domain 
wall fermions]^. However, calculations using im- 
proved Wilson fermions were also reported at this 
meeting and it is still unclear which method 
will prove most advantageous. 
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In this study, we present preliminary results for 
the K ^ n and X — > matrix elements of these 
operators, which when combined with lowest or- 
der chiral perturbation theory, yield the desired 
K TTTT matrix elements Q. The contribution of 
R. Mawhinney in these proceedings takes up this 
point Here we are concerned with the simpler 
K ^ TT and K ^ lattice matrix elements. Also, 
the CP-PACS collaboration has presented a very 
similar calculation at this meeting]^]. 

2. Theoretical Framework 

The AS* = 1 effective weak Hamiltonian is gen- 
erated from the fundamental Standard Model La- 
grangian by integrating out the top quark and 
W boson. The resulting effective Hamiltonian is 
then evolved to a low scale (/i <ti Mw) appropri- 
ate for lattice calculations using the renormaliza- 
tion group equations. The effective Hamiltonian 
above the charm threshold is: 
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where fi is the renormalization scale, Q are a 
basis of local four-quark operators which are 
closed under renormalization, C'i{ii) the cor- 
responding Wilson coefficients, and Vqq' the 
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Cabibbo-Kobayashi-Maskawa mixing matrix ele- 
ments which are fundamental parameters of the 
Standard Model. 

The effective operators renormalized at the 
scale /i are 

Ql = SaluPLdaUp-fyPLUp 

= SalvPhda ^ qpltiliyP(L,R)q(3 
u,d,s,c... 

= SalvPhdp ^ qi3ltJ.lvP(L,B.)qa 
u.d,s^c... 
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Qlc = SalvPhdaCp^^PLCfj 
Q2c = SaluPLdpCp-fuPhCa, 
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with color indices a and /3, P(l,r) = 1 T 75j 
the sums are taken over active quark flavors at 
the scale /x, and summation over v is implied. 
Ql, 2,1c, 2c are often referred to as current-current 
operators, Q3-6 QCD penguin operators, and 
Qt-iq electroweak penguin operators. It is use- 
ful to split the above operators according to their 
isospin, Qr = Qf'''^ +Qf'^\ 

As mentioned earlier, lowest order chiral per- 
turbation theory relates K tttt matrix elements 
to a linear combination of K it and X — s- 0. 
For all Qi except the electroweak penguins Qj^s 
and all particles at rest, we haveQ 
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where niM is the meson mass for unphysical pseu- 
doscalar states with vris = md- For Qi which 
transform in a (27,1) chiral multiplet a2 — 0. 
Note that each matrix element vanishes linearly 
with the meson mass squared. This is an im- 
portant prediction of chiral perturbation theory. 



and therefore QCD, and provides a solid test of 
the chiral symmetry properties of domain wall 
fermions. The strength of the above approach 
is that it allows less computationally demanding 
K IT and A' matrix elements to be cal- 
culated on the lattice. A significant drawback to 
this approach is that it manifestly does not con- 
tain information on the final state interactions of 
the pions (for calculation of e' the final state s- 
wave scattering phases from experiment can be 
put in by hand, however). 

In the case of the electroweak penguins Qj^s^ 
the contribution in lowest order chiral perturba- 
tion theory to the K ^ tt matrix elements is con- 
stant in the chiral limit. The K ^ matrix ele- 
ment, however, still vanishes. 

Since on the lattice 02 is quadratically diver- 
gent, the process of combining the second and 
third lines in Eq. 4 to obtain an expression for 
K ^ TTTT requires a delicate cancellation of this 
quadratic divergence. See R. Mawhinney's con- 
tribution for details. 

Finally, lattice counterparts of the operators in 
Eq. 1^ must be matched to the continuum and 
renormalized since they are logarithmically di- 
vergent after power divergences have been sub- 
tracted. In addition, operators in the same 
symmetry multiplets mix through renormaliza- 
tion group running from M\y down to the low 
scale /i. This poses a serious challenge for lat- 
tice calculations. In our calculation this prob- 
lem is handled remarkably well with the nonper- 
turbative renormalization method of the Rome- 
Southampton group which was explained in the 
talk by C. Dawsong. 

3. Simulation details 

We have calculated matrix elements on 200 
quenched gauge configurations at /3 = 6.0, with 
lattice four volume 16^ x 32, domain wall fermion 
extra dimension size Lg = 16, and domain wall 
height Afs = 1.8. We have calculated with light 
quark masses m/ = 0.01 — 0.05, and charm 
quark masses rric = 0.1 — 0.4. The physical kaon 
state made from degenerate quarks corresponds 
to nif Ri 0.02, and nic ~ 0.5 for the physical 
charm quark. 
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We extract matrix elements from three-point 
correlation functions. The external pseudoscalar 
states are interpolated from wall sources near the 
time direction boundaries, t = 5 and 27, and 
the operator is inserted between them. When 
the operator is far from either boundary, the de- 
sired lowest mass states dominate the correlation 
function. The forward and backward (in time) 
quark propagators used to interpolate the K and 
TT states are linear combinations of propagators 
computed with periodic and anti-periodic bound- 
ary conditions which amounts to doubling the 
gauge field configuration in the time direction. 
The closed fermion loops necessary for operators 
that have self contractions are computed from 
a complex Gaussian random source spread over 
time slices 14-17. All results are given as aver- 
ages over these four time slices. 

We have performed several important checks of 
our computer code. Most importantly, a com- 
pletely independent check code was written to 
compare with our two production versions (the 
check code and one of the production codes run 
on the QCDSP supercomputer and the other pro- 
duction code, based on the MILC code, runs on 
the NERSC T3E). Output from each code gen- 
erated on the same configuration agreed up to 
machine precision. 

As a final useful check, the left-left operators 
in Eq. |^ go into themselves under a Fierz trans- 
formation. Thus color-mixed contractions can be 
compared to corresponding color-diagonal ones. 
We find perfect agreement in all cases. 

The following results were obtained on the 
RIKEN BNL and Columbia University QCDSP 
supercomputers. 

4. Results 

Fig. I shows (7r|Q^^/^^|X) as a function of 
quark mass, w/ = 77is — m^. An uncorrelated 
linear extrapolation yields a zero intercept, within 
statistical errors, which is in agreement with chi- 
ral perturbation theory. For strictly low energy 
QCD observables, we expect quantities to vanish 
ai irif = —rrires^- Since the AI — 1/2 opera- 
tors have contributions from physics scales near 
the (high energy) lattice cut-off, this is no longer 



true. Thus, for this matrix element, the statisti- 
cal errors are not small enough to resolve these 
systematic effects. However, in R. Mawhinney's 
contribution, we see that such effects are visible 
in the subtracted operator^. Presumably this is 
due to the strong statistical correlations between 
the A' — > TT matrix element and the subtraction 
term. In Fig. || we show a similar plot for Qq. 
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Figure 1. The K it matrix element of the bare 
operator Q^2^'^^ . 



Here, an uncorrelated linear extrapolation has a 
non-zero intercept of roughly three standard de- 
viations. Note that it vanishes for m/ > 0. Thus 
explicit chiral symmetry breaking effects are visi- 
ble, though small. In Fig. ||we show (ttIQ^/^^IA-) 
which exhibits noticeable nonlinearity and does 
not vanish in the chiral limit, as expected. 

Finally, we show an example of a X — > ma- 
trix element in Fig. |. Note that for the AT ^ 
matrix elements we use ^ and fit to the 
form i^\Qi\K) I (^\s^^d\K) = const + {rus — md)r]i 
where chiral perturbation theory predicts const — 
0. This ratio is useful since it is exactly the coef- 
ficient of the subtraction operator used to remove 
the quadratic divergence in {Tr\Qi\K). Since the 
quark masses enter as a difference, we expect ex- 
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Figure 2. The K ^ tt matrix element of the bare Figure 3. The K ^ tt matrix element of the bare 

operator Qq (charm contribution not included) . operator Q^g^^^ . 



plicit chiral symmetry breaking effects which do 
not depend on the quark mass to cancel. From 
the fit depicted in Fig. ^ wc find that the constant 
term is zero within errors. Note that this ratio is 
extremely well resolved, and quite linear. 

We take this opportunity to quote our value for 

the kaon B parameter, B'j^'^\2 GeV) = 0.538(8), 
with Zll/Z\ = 0.928(6) computed nonperturba- 
tively in the RI scheme and matched to the MS- 
NDR scheme 0. The errors are statistical only, 
and the error on Bk is obtained by adding the er- 
rors on the matrix element and the renormaliza- 
tion factor in quadrature. Our value is lower than 
the one quoted by Taniguchi at this meeting||], 
probably due to the fact that our renormalization 
constant is lower than the pcrturbative one used 
in that study. Note our result is for f3 — 6.0 with 
a single lattice size. 
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Figure 4. The K ^ matrix element of the bare 
operator Qq (charm contribution not included). 
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